Mon. Not. R. Astron. Soc. 0Q0,[Tiri31(2Q02) Printed 8 January 2010 (MN BOgX style file v2.2) 



Near-IR dust and line emission from the central region of Mrk 1066: 
Constraints from Gemini NIFS 

Rogemar. A. Riffel 1 ' 2 *, Thaisa Storchi-Bergmann 2 and Neil M. Nagar 3 

1 Universidade Federal de Santa Maria, Departamento de Fi'sica, Centro de Ciencias Naturals e Exatas, 97105-900, Santa Maria, RS, Brazil 
2 Universidade Federal do Rio Grande do Sul, Departamento de Astronomia, Institute de Fi'sica, CP 15051, 91501-970, Porto Alegre, RS, Brazil 
3 Astronomy Group, Departamento de Fi'sica, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 



Accepted 1988 December 15. Received 1988 December 14; in original form 1988 October 11 



ABSTRACT 

We present integral field spectroscopy of the inner 700 x 700 pc 2 of the Seyfert galaxy 
Mrk 1066 obtained with Gemini's Near-Infrared Integral Field Spectrograph (NIFS) at a spa- 
tial resolution of *35 pc. This high spatial resolution allowed us to observe, for the first time 
in this galaxy, an unresolved dust concentration with mass ~ 1.4 x 10 -2 M Q . This unresolved 
concentration, with emission well reproduced by a blackbody with temperature ~830K, is 
possibly part of the nuclear dusty torus. We compared maps of emission-line flux distributions 
and ratios with a 3.6 cm radio-continuum image and [Om] image in order to investigate the 
origin of the near-infrared emission. The emission-line fluxes are elongated in PA= 135/315° 
in agreement with the [O m] and radio images and, except for the H lines, are brighter to the 
north-west than to the south-east. This close association with the radio hot spot implies that at 
least part of the emitting gas is co-spatial with the radio outflow. The H emission is stronger 
to the south-east, where we find a large region of star-formation. The strong correlation be- 
tween the radio emission and the highest emission-line fluxes indicates that the radio jet plays 
a fundamental role at these intensity levels. At lower emission-line fluxes this correlation dis- 
appears suggesting a contribution from the plane of the galaxy to the observed emission. The 
H2 flux is more uniformly distributed and has an excitation temperature of = 2100K. Its ori- 
gin appears to be circum-nuclear gas heated by X-rays from the central active nucleus. The 
[Fe 11] emission also is consistent with X-ray heating, but its spatial correlation with the ra- 
dio jet and [O m] emission indicates additional emission due to excitation and/or abundance 
changes caused by shocks in the radio jet. The coronal-line emission of [Ca vm] and [S ix] are 
unresolved by our observations indicating a distribution within 18pc from the nucleus. The 
reddening map obtained via the Pa/3/Bry line ratio ranges from E(B - V) ~ to E(B - V) ~ 1 .7 
with the highest values defining a S-shaped structure along PAa:135/315 . The emission-line 
ratios are Seyfert-like within the ionization cone indicating that the line emission is powered 
by the central active nucleus in these locations. Low ionization regions are observed away 
from the ionization cone, and may be powered by the diffuse radiation field which filters 
through the ionization cone walls. Two regions at 075 south-east and at 1" north-west of the 
nucleus show starburst-like line ratios, co-spatial with an enhancement in the emission of the 
H lines. We attribute this change to additional emission from star forming regions. The mass 
of ionized gas is M HII sl.7x 10 7 M G and that of hot molecular gas is M Hl =3 3.3 x 10 3 M o . 

Key words: galaxies: individual (Mrk 1066) - galaxies: Seyfert - galaxies: ISM - infrared: 
galaxies 



1 INTRODUCTION 

The study of the gas distribution and excitation in the Narrow Line 
Region (NLR) of active galaxies has a fundamental importance in 
the understanding of the physical phenomena in the vicinity of 

* E-mail: rogemar@smail.ufsm.br 



the super-massive black holes (SMBH) in their centers. The ex- 
citation of the inner NLR can reveal how the radiation and mass 
outflows from the active galactic nucleus (AGN) interact with the 
surrounding inter-stellar medium (ISM). Several studies have been 
aimed to investigate the line emission from the NLR gas (e.g. 
Rodriguez- Ardila et al. 2004; Rodriguez-Ardila, Riffel & Pastoriza 
2005al : iRiffel et alj l200fj, |2008l , l2009al ; IStoklasova et all 120091; 
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Figure 1. Top left: Large scale image of Mrk 1066 obtained at the R-band with the 1.22 m Oschin Telescope on Palomar Mountain jLasker et alj 1996). The 
central box shows NIFS field of view. Top right: 2.22//m continuum map from NIFS observations. Bottom: Spectra for the nuclear position (N at top-right 
panel) and for the position A ((X'5 north-west from the nucleus) for the J and Kj-bands for an aperture of (X'25x'.'25. The emission lines are identified at the 
nuclear spectrum and in the extra-nuclear spectrum we mark the K-band CO absorption band heads. 



IStorchi-Bergmann et al1l2009l ; iBarbosa et alj|2009l) . Nevertheless, 
its physics is still not completely understood. 

Besides the line emission, the nature of the nuclear continuum 
emission is an important issue in the study of AGNs. In particular, 
the study of near-infrared (hereafter near-IR) continuum can be 
used to probe the emission from hot d ust located in the putative 
torus postulated by the Unified Model dAntonuccill 19931) . Several 
studies revealed the presence of hot dust with temperatures of 
T ~ 900 - 1500 K and m as ses of M HD ~ 10~ 6 - 10~ 2 M, (e.g 



Rieke & Lebofsky 1981; Rodriguez- Ardila, Contini & Viegas 



2005b): iRodrfguez-Ardila & Mazzalavl l2006t iRiffel et al 
2009a b e). The high spatial resolution provided by adaptive 
optics spectroscopy in large telescopes can be used to constra int 
the location of the near-IR emitting dust (e.g. lRiffel et al .l2009bh . 

This work is aimed to investigate the near-IR emission- 
line flux distributions, gas excitation and continuum emission 
within the inner kiloparsec of Mrk 1066 from spectroscopy ob- 
tained with t he Gemini's Near-Infra red Integral Field Spectro- 
graph (NIFS, iMcGreeor etatl [2003). In a following paper we 
discuss the kinematics of both the gas and stars. Mrk 1066 was 



selected for this study because it presen ts extended radio (e.g 
lUlvestad & Wilson 1989; Nagar et al. 1999) an d near-IR line emis- 
sion ( e.g. IRodrfguez-Ardila. Riffel & Pastor iza 2005a; iKnop et al] 
2001) allowing us to investigate the role of the radio jet in the emis- 
sion of the NLR gas. In addition, its central kp c has a high extinc- 
tion in the optical (e.g. IStoklasova et alll2009l) and thus a near-IR 
study is better suited for the investigation of the gas excitation and 
kinematics. 

Mrk 1066 is an SBO+ galaxy harboring a Seyfert 2 nucleus; it 
is located at a distance of 48.6 MpcQ for which 1" corresponds to 
235 pc at the galaxy. Hubble Space Telescope (HST) narrow-band 
images show a "jetlike" feature in the [Om]+Uy6 emission extend- 
ing up to 1'.'4 north-west from the nucleus along the position an- 
gle PA=315°, w hile the Hq+[Nn] image is extended to both sides 
of the nucleus dBower et al] I1995T) . Radio continuum images of 
Mrk 1066, at 3.6, 6 and 20 cm, show extended emission up to 1'.'5 to 
both sides of the nucleus oriented approximately along the same PA 



1 Distance quoted in NASA/IPAC Extragalactic Database (NED 
http://nedwww.ipac.caltech.edu I 
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Figur e 2. Continuum images o f the central 3"x3" of Mrk 1066. Top left: Optical continuum image obtained with the HST-WFPC2 using the filter F606W 
I Malka n. Goriian & TanJ 19981) . Top right: H-band continuum image obtained with HST-NICMOS with filter F160W. Bottom: Continuum images centered at 
1 .23 fjm (left) and 2.22 /jm (right) from NIFS observations. The HST images are in arbitrary units and NIFS images are in logarithmic flux units (erg s _I cm -2 
A ). The straight line shows the orientation of the major axis of the galaxy and the thick contours overlaid to the K-band map are from the optical continuum 
image. 



of the optical line emission iUlvestad & WilsonH l989; Nag ar et all 
1 19991) . Long slit spectroscopy shows that the near-IR emission lines 
are extended up to 5" from the nucleus along the PA= 1 35/3 15° with 
distinct flux distributions for [Fen], H 2 and H recombinat ion lines 
suggesting different emission processes tenopetal]200lh . 

This paper is organized as follows: In Section 2 we present 
the description of the observations and data reduction. In Sec. 3 we 
present the results for the continuum and line emission, as well as 
line-ratio maps. The discussion of the results is presented in Sec. 4 
and the conclusions are presented in Sec. 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The Gemini NIFS data 

The Integral Field Unit (IFU ) spectroscopic data of Mrk 1066 were 
obtained with Gemini NIFS (McGregor et al. 2003) operating with 
the Gemini North Adaptive Optics system ALTAIR in Septem- 
ber 2008 under the programme GN-2008B-Q-30. The NIFS has 



a square field of view of as 3'.'0 x 3'.'0, divided into 29 slices with 
an angular sampling of 0'.' 103x0'.' 042, optimized for use with AL- 
TAIR. 

The observing procedure followed the standard Object-Sky- 
Sky-Object dither sequence, with off-source sky positions since the 
target is extended, and individual exposure times of 600 s. Two set 
of observations were obtained at different spectral ranges: the first 
at the J-band, centered at 1.25 fim and covering the spectral region 
from 1.15yum to 1.36yum, and the second at the Kj-band, centered 
at 2.3 yum covering the spectral range from 2. lOyum to 2.26 yum. At 
the J-band, the selected instrument configuration was the LG5603 
grating and ZLG0601 filter resulting in a spectral resolution of ~ 
1.7 A, as obtained from the measurement of the full width at half 
maximum (FWHM) of arc lamp lines. The Kpband observations 
were done using the KLG5607 grating and HK.G0603 filter and 
resulted in a spectral resolution of FWHM» 3.3 A. 

The total exposure time at each band was 4800 s, consisting of 
8 individual on-source exposures. A small spatial dithering of 071 
was applied between on-source exposures in order to correct the 
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data from bad pixels. The FWHM of the spatial profile of the star 
was a' 13±O702 for the J-band and a' 15±a'03 for the K r band. 

The data reduction was accomplished using tasks contained in 
the nifs package which is part of gemini iraf package as well as 
generic iraf tasks. The data reduction followed the standard pro- 
cedure, which includes the trimming of the images, flat-fielding, 
cosmic ray rejection, sky subtraction, wavelength and s-distortion 
calibrations. In order to remove telluric absorptions from the galaxy 
spectra we observed the standard star HIP 15925 just after the J- 
band observations of the galaxy and HIP 10559 just before the K r 
band observations. The galaxy spectra was divided by the normal- 
ized spectrum of the tellurc standard star using the nftelluric task 
of the nifs.gemini.iraf package. The galaxy spectra were flux cal- 
ibrated by interpolating a blackbody function to the spectrum of 
the telluric standard and the J and K -bands data cubes were con- 
structed with an angular sampling of 0705 x0705 for each individual 
exposure. Finally the individual data cubes were combined using a 
sigma clipping algorithm in order to eliminate bad pixels and re- 
mainding cosmic rays by mosaicing the dithered spatial positions. 
The final J and K data cubes contain about 4200 individual spectra 
and cover the central 3"x3", corresponding to ^700x700 square 
parsecs at the galaxy. 

2.2 The 3.6 cm radio image 

The 3.6 cm Very Large Array (VLA) radio image was made using 
data ta ken on 31 December 1992, initial ly published inlNagar et al.l 
( 1999), and reprocessed as described in Mund ell et all | |2009|) . 



3 RESULTS 

In the top-left panel of FigJj] we present a large scale R-bmd 
image of Mrk 10 66 from the Palomar Observatory Sky Survey 
(Laske r1-t~aill996h . In the top-right panel we present an image ob- 
tained from the NIFS data cube for the continuum around 2.22 /im. 
In the bottom panels we present two characteristic IFU spectra: 
the nuclear one and a spectrum from a location at 075 north- 
west of the nucleus. In the J-band we identified the [P n] emis- 
sion lines at 1.14713 and 1.18861 fan, the [Fen] emission lines at 
1.25702, 1.27912, 1.29462, 1.29812, 1.32092 and 1.32814pm, the 
Hi Pa/3 A 1.282 16pm, the Hen line at 1.16296 /mi and the [Six] 
coronal line at 1.25235 fim at the nucleus. The H2 emission lines at 
2.12183, 2.15420, 2.22344, 2.24776, 2.40847, 2.41367, 2.42180, 
2.43697 and 2.45485 fim are identified in the K-band spectra, as 
well as the Hi Bry /i2.16612/mi, the He iA2. 14999 fim and the 
[Ca vm] A2. 32204 coronal line at the nucleus. In the K-band 
spectra we have also identified the CO stellar absorption band heads 
around 2.3 fim. 

3.1 The continuum emission 

Figure|2] presents optical and near-IR continuum images of the 
central 3"x3" of Mrk 1066. The top left panel shows an opti- 
cal image obtained with the HST Wid e Field Planetary Cam- 
era 2 (WFPC2) through the filter F606W jMalkan. Goriian & Tard 
1998). This image presents two spiral arms oriented approximately 
along PA=128° (approximately coincident with the major axis of 
the galaxy), which extend up to 172 to both sides of the nucleus. 
Along the spiral arms there are several knots with enhanced emis- 
sion. This image is i n good agreement with the optical continuum 
images presented bv lBower et alj dl995h . 
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Figure 3. K/J continuum ratio obtained by averaging continuum regions 
(in units of erg s cm~ 2 A - ' ) centered on 2.22 fim and 1 .23 fim. The thick 
contours are from the reddening map for Fig. [5] 

In the top right panel of Fig.[2] we present an HST H-band 
image obtained with the Near Infrared Camera and Multi-Object 
Spectrometer (NICMOS) through the filter F160W (proposal ID 
7330 by J. S. Mulchaey).In the bottom panels of Fig.[2]we present 
the J and K-band continuum images obtained from the NIFS data 
cubes by averaging continuum regions centered at 1.23 yum and 
2.22yum.The near-IR continuum is more extended along PA=128°, 
which is similar to the orientation of the spiral arms observed in the 
optical continuum and approximately the orientation of the major 
axis of the galaxy. The spiral arms can also be seen in the near-IR, 
although with "less contrast" than in the optical, probably due to 
the stronger effect of the dust and contribution of line emission in 
the optical. A comparison between the J and K continuum shows 
that the extra-nuclear emission is bluer than the nuclear one, indi- 
cating the presence of dust nucleus. This behavior can be observed 
in Fig. [3] which presents a K/J continuum ratio map obtained from 
the NIFS data cubes by averaging the flux within spectral window 
of 50 A in regions free of emission and absorption lines centered at 
2.22 pm and 1.23/um. 

3.2 Emission-line flux distributions 

As shown in Fig. Q] several emission lines can be observed in the 
spectra. Table[T]presents the measured fluxes for the emission lines 
within an aperture of 0725 x0725 at the four positions identified in 
the top-right panel of Fig.Q] the nucleus, position A at 075 north- 
west from the nucleus, position B at 075 south-east and position C 
at 075 north-east from the nucleus. Most emission lines peak at the 
nucleus, except for [Fe 11] and [P 11] which peak at 07 5 north-west of 
the nucleus. 

In order to map the gas emission we have integrated the flux 
under the emission-line profiles and subtracted the underlying con- 
tinuum obtained from two spectral windows, one from each side 
of the profile. Figure|4] presents the resulting flux maps. For most 
emission lines, the highest intensity levels most extended along 
PA» 135/315°, which is the orientation of the radio structure. 
A detailed inspection of each panel reveals distinct structures for 
each ionization level. The [Pn], [Fen] and H2 emission lines are 
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Figure 4. Integrated flux maps for the [Pn]/U.8861 ^m, [Fen]/U.2570/im, Pa/3, [Fen]^1.3109/an, H 2 /12.1218,um, Bry, [Cavm]i2.3220^m and 
H2/i2.1218/im emission lines in units of 10~ 17 ergs -1 cm -2 . The central cross marks the position of the continuum peak. The thick contours overlaid on 
the Pa/3, [Fen]/U.2570/mi and H2/12.1218/im maps are from the 3.6 cm radio continuum image. 
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Table 1. Measured emission-line fluxes for the four positions marked in Fig.[T]witmn (X'25x(X'25 aperture in units of 1CT 16 erg s -1 cm 



A vac (A) 


ID 




Nucleus 


A 


B 


C 


11471.3 


[Pii]'D 3 - 3 Pi 




5.07±2.23 


2.75±0.71 


1.50±0.47 




11629.6 


Hen7-5 




2.41±1.52 


0.99±0.32 


0.67 ±0.13 




11886.1 


[Pn] 'D 2 - 3 Pi 




6.21±0.69 


6.54±0.87 


2.80±0.38 


0.41±0.28 


12523.5 


[Six] 3 /-, - 3 P 2 


0.60±0.37 








12570.2 


[Fen] - 


a 6 D 9n 


18.88±0.80 


21.15±1.82 


13.06±0.95 


5.76+1.21 


12706.9 


[Fen] a 4 Di/2 — 


a 6 D 1/2 


0.52±0.41 


n Qfi+o 38 


0.43±0.17 


0.14+0. 1 1 


12712.1 


[Feiv] 2 G^ /2 - 2 


D 5 "> 


0.21+0.19 








12791.2 


[Fen] a 4 Z>3/2 - 


a 6 D i/2 


0.94±0.35 


1.34±0.65 


1.27±0.32 


0.36±0.25 


12821.6 


HiPa/3 




41.18±4.26 


25.17+1.98 


26.41±1.98 


6.79±0.55 


12946.2 


[Fen] a 4 Ds/2 - 


a 6 D 5l2 


1.91±0.87 


1.71+0.16 


0.82±0.29 


0.70±0.13 


12981.2 


[Fen]a 4 Z) 3/2 - 


a 6 D ul 


1.72+0.55 


0.53±0.22 


0.94±0.23 


0.46±0.29 


13209.2 


[Fen] a A Dj/2 - 


a 6 D 1/2 


4.99±0.73 


6.32±0.64 


4.03±0.52 


2.21±0.26 


13281.4 


[Fen]u 4 Z) 5 /2 - 


a b D vl 


1.54±1.05 


1.31+0.59 


0.54±0.40 




21218.3 


H 2 2 1-0S(1) 




13.41±1.02 


11.74±0.49 


5.27±0.48 


2.52+0.11 


21499.9 


He 1 3 S i - 3 P 1 ! 




0.40±0.21 


0.11 ±0.21 


0.27±0.19 


0.18±0.07 


21542.0 


H 2 1-0 S(2) 




0.51±0.23 


0.45±0.50 


0.25 ±0.18 


0.16±0.04 


21661.2 


Hi Bry 




17.14±1.16 


10.59±0.60 


9.81±0.53 


1.28±0.08 


22234.4 


H 2 1-0 S(0) 




2.99±0.92 


2.63±0.64 


1.96±0.72 


0.65±0.06 


22477.6 


H 2 2-1S(1) 




1.47±0.62 


1.17+0.38 


0.67±0.33 


0.33±0.02 


23220.4 


[Cavm] 2 P° /2 - 


2 pO 

M/2 


4.64±1.42 








24084.7 


H 2 1-0 Q(l) 




10.04±1.33 


8.14±0.45 


4.82±0.86 


1.42+0.11 


24136.7 


H 2 1-0 Q(2) 




3.03±0.47 


2.47±0.42 


1.47±0.41 


0.49±0.10 


24218.0 


H 2 1-0Q(3) 




10.61±1.05 


8.32±0.57 


5.23±0.56 


1.50±0.21 


24369.7 


H 2 1-0 Q(4) 




2.73±1.45 


2.26±0.38 


1.30±0.48 


0.47±0.09 


24548.5 


H 2 1-0 Q(5) 




7.08±2.18 


5.41±0.43 


2.82±0.77 


4.36±2.24 



~2 times brighter to the north-west than to the south-east, while 
Pa/3 and Bry are approximately equally bright to both sides of the 
nucleus. The peak fluxes of the [Pn] and [Fen] emission lines 
occur at «(y.'85 north-west of the nucleus, approximately at the 
same position where there is a hot spot in the radio map. The H i 
recombination-line emission peaks at the nucleus and presents a 
secondary peak at xQf!5 south-east from it. Another difference be- 
tween the Hi and [Fen] emission lines is that the flux distribution 
of the latter traces the radio structure better than the former. The 
H 2 emission also presents a peak at the nucleus and a secondary 
peak at aC.' 65 north-west of the nucleus, in a region between the 
nucleus and the radio hot spot. The emission in the coronal lines of 
[Cavm]/12.3220yum and [Six] Al. 2524 /urn are unresolved by our 
observations. 



3.3 Line-ratio maps 

In order to investigate the excitation mechanisms of the [Fe n] and 
H 2 emission lines we constructed flux ratio maps, shown in Fig- 
ure[5] In this figure we also present a reddening map obtained from 
the Pa/3/Bry line ratio as 



E(B -V) = 4.74 log 



where F P 



5.88 \ 

F p a p /F Bry J 



(1) 



' PaB and F Bly are the fluxes of Pa/3 and Bry emis- 
sion lines, respectively. We hav e used the reddening law of 
ICardelli. Clayton & Mathisl Jl989h and adopted the intrinsic ra- 
tio F PaS /F Brv = 5.88 cor responding to case B recombination 
dOsterbrock & FerlandkOOrjh . The resulting E(B - V) map is shown 
in the top left panel of Fig.[5] The highest values are observed in 
a S-shaped structure, resembling the spiral arms observed in the 
continuum image overlaid as thick contours on the reddening map. 



Typical values within this structure are E(B-V) ss 1.0, with several 
knots of higher values reaching E(B - V) ~ 1.7. In regions away 
from the center of the S-shaped structure the reddening decreases to 
E(B - V) a 0.4 or bellow. The errors in E(B — V) are approximately 
0.15. 

A comparison between the emission-line reddening and the 
continuum ratio map is shown in Fig. [3] The thick contours, rep- 
resenting the reddening map, show that the S-shaped structure is 
also present in the continuum ratio map as a region of higher K/J 
ratio. In order to estimate the reddening in the continuum, we note 
that, from Fig. [3] a continuum ratio (K/J)i »0.4 corresponds to the 
locations with E{B — V)\ » 0.4 in the gas. If the higher (K/J) 2 ratio 
along the S-shaped structure is due to reddening, the corresponding 
E(B - V) 2 can be calculated from 



X [E(B - Vh - E(B - V) 2 ] = log 



(K/J), 



(K/J) 2 



(2) 



where x = -0A(A K /A v - Aj/A V )R V = 0.2083, where (A K /A V ) = 
0.114, (Aj/A v ) = 0.282 and R v = 3.1 for the extinction law of 
ICardelli. Clayton & Mathisl dl989l) . Using (K/J) 2 = 0.5 measured 
in the S-shaped structure (from Fig. [3} we obtain reddening for 
the continuum of E(B - V) 2 ~ 0.9, which is somewhat smaller 
than the one observed for the emitting gas (top-left panel of Fig.|5]l. 
Nevertheless, right at the nucleus the K/J ratio is much higher than 
in the rest of the S-shaped structure, what indicates that it is not due 
to reddening. The nature of this redder continuum will be discussed 
in SecRTTI 

In the top-right panel of Fig.[5] we present the 
[Feii]/ll.2570//m/Pa/3 intensity ratio, which can be used 
to investigate the excitation mechanism of the [Fen] 
emission. Seyfert galaxies present ty pical value s for 
this ratio between 0.6 and 2.0. (e.g. iRiffel et alj l2006t 
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Figure 5. Top: Reddening map obtained from the Pa/?/Bry line ratio (left) and [Fen]/11.2570^m/Pa/3 line ratio map (right). Bottom: 
[Fen]/ll. 2570 yum/[Pn]/ll. 8861 fim line-ratio map (left) and H2/12.1218 J um/Bry ratio (right). The thick contours overlaid to the [Fen]/Pa/3 and H2/Bry maps 
are from the radio image. 



Rodriguez- Ardila et al. 2004; Rodrigu ez- Ardila. Riffel & Pastoriza] 



2005a; IStorchi-Bergmann et alj l2009f) . The smallest value ob- 
served in Mrkl066 is [Fen]^1.2570^m/Pa/? = 0.33 ± 0.03 at 
073 south-east from the nucleus, while the highest values are 
observed in regions located at distances larger than 1", reaching 
values of up to 1.5. 

Another line ratio that can be used to investigate the [Fen] 
excitation mechanism is [Fe n]/tl. 2570 //m/[Pn ]/ll. 8861 /an (e.g. 
Oliva et al ] |200ll; IStorchi -Bergmann et al . 2009). We present this 
ratio in the bottom-left panel of Fig.[5] The lowest values of w 3 
are observed at the nucleus and in regions to north-west at small 
distances from the nucleus, while the highest values up to 9.5 are 
observed at larger distances from the nucleus. 

In the bottom-right panel of Fig.[5] we present the 
H2/l2.1218/jm/Bry ratio map, which is useful to investigate the 
excitation of the H 2 emitting gas. This ratio pre sents values rang- 
ing from 0.6 to 2 for Seyfert galaxies (e.g. |Riffel et al.l [2006]; 
Rodriguez- Ardila et al. 2004; Rodriguez- Ardila, Riffel & Pastoriza 
2005al ; IStorchi-Bergmann et alj|2009h . In Mrk 1066, the lowest 
values of 0.29 ± 0.06 are observed at 1" north-west from the nu- 
cleus and the highest values of up to 2.6 are observed at distances 



of >077 in regions away from the radio structure. Another region of 
small values (« 0.4) is observed at »075 south-east of the nucleus, 
approximately coincident with the secondary peak observed in the 
H i recombination-lines flux maps (see Fig.[4), while at the nucleus 
H 2 A2.1218yum/Br7 = 0.75 +0.07. 



4 DISCUSSION 
4.1 Dust emission 

As observed in Fig.|3](see also Figs.[TJand[2) the extra-nuclear con- 
tinuum is bluer than the nuclear continuum. At the nucleus the 
ratio between the continuum at 2.22 jim and at 1.23 yum is « 1, 
while at larger distances from the nucleus, typical values are « 0.45 
along PAr; 128° and even smaller at locations away from this PA 
(< 0.38). In order to to investigate the nature of the K-band excess 
at the nucleus, we have extracted a nuclear spectrum within an aper- 
ture of 0725 radius, which includes the nuclear flux down to 20% 
of the peak continuum flux in the K band. An extra-nuclear spec- 
trum for a circular ring within radii 0725< r <0735 was extracted 
to represent the stellar population. This spectrum was normalized 
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according to the ratio of the nuclear/extra-nuclear extraction aper- 
tures. Both spectra were corrected for extinction using the redden- 
ing values obtained in Sec. 13.31 E(B — V) » 1.03 for the nucleus 
and E(B — V) ~ 0.69 for the stellar population. 

Under the assumption that the underlying stellar population in 
the nuclear spectrum is the same as in the extra-nuclear spectrum, 
we have subtracted the latter from the former. The resulting nuclear 
continuum was then modeled by the sum of two components: a 
blackbody function, which dominates the emission in the K band 
and a power-law, which is important for the J-band emission. The 
fit was performed using the nfitld task of the stsdas iraf package 
selecting only spectral regions not affected by line emission. The 
best fit was obtained for a power-law F A oc A~ 2 A±0 3 and a blackbody 
function with temperature T = 863 ± 30 K. 

We then considered the possibility that the contribution of the 
stellar population to the nuclear spectrum is higher than the one 
derived only from the normalization of the fluxes to the same ex- 
traction aperture. This possibility is suggested by the fact that the 
stellar CO absorption bands in the nuclear spectrum do not disap- 
pear after the subtraction of the extra-nuclear spectrum. In order to 
completely eliminate the stellar absorptions the extra-nuclear spec- 
trum had to be multiplied by a factor of 1.45. After this correction 
we obtain F A cc A- 4J±0A and T = 829 ± 30 K. In the top panel of 
Fig. [6] we show the nuclear and extra-nuclear spectra corrected by 
extinction with the latter properly normalized as described above. 
The bottom panel shows the subtracted spectrum, as well as the re- 
sulting fit as a continuous line. The blackbody function is shown as 
a dashed line a nd the powe r -law as a dotted line. 

Recently, iRiffel et al.l |2009cJl presented stellar population 
synthesis for a sample of Seyfert galaxies, including Mrk 1066, us- 
ing near-IR spectra extracted within an aperture of 1'.' 6x078. Be- 
sides synthetic Simple Stellar Population (SSP) models, the au- 
thors considered the contribution of a power-law to represent the 
featureless continuum from the AGN and five blackbody func- 
tions for temperatures ranging from 800 K to 1400 K in order 
to account for possible dust emission. They concluded that the 
dust emission is not necessary in the case of Mrk 1066. Never- 
theless, as shown in Fig.|6l a K-band exce ss is clearly present in 
our nuclear spectrum. IRiffel et all |2009c) did not identified this 
component probably due to their larger aperture, which covers 
an area more than 6 times ours. Unresolved dust emission from 
the nucleus of Seyfert galaxies is commonly observed and is usu- 
ally attributed to the dusty torus postulated by the Unified Model 
(e.g. IRiffel et all l2009allbld; iRodrfguez-Ardila & Mazzalavl 120061 : 
Rodriguez-Ardila, Contini & Viegasll2005bl) . Previous works have 
shown that the dust emission in the near-IR for Seyfert 1 galax- 
ies give blackbody temperatures T > 1200 K, while for Seyfert 2 
galaxies the blackbody temperature is usually T < 900 K, in agree- 
ment with our result for Mrk 1066. This difference can be inter- 
preted as due to the fact that in Seyfert 1 galaxies one can see 
the inner wall of the torus, where the temperatures are c lose to the 
evaporation temperature o f graphite and silicate grains dBarvainisl 
Il987l; iGranato & Daneselll994l) . In Seyfert 2 galaxies we do not 
see the inner wall, only dust located at larger distances from the 
nucleus, heated by a more diluted radiation field, part of which is 
also blocked by clouds in the inner region of the torus (e.g. Elitzur 
2008). Since the K-band continuum of Mrk 1066 is unresolved by 
our observations, we can constraint the location of the emitting dust 
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to with in the inner w 18 pc, although, as discussed in Rif fel et al.l 
(2009b), the dust is probably located much closer in. 

We can derive the mass of the emitting dust from the parame- 
ters obtained from the fit of the blackbody function, together with 
assumptions about the physical properties of the dust grains. We es- 
timate the infrared luminosity of each grain assuming that the dust 
composition is graphite, by 

Z/; ir = 4;rVe v B v (7 gr ) [ergs-'Hz- 1 ] (3) 

where a is the grain radius, Q y = q ir V is the absorption effi- 
ciency and S r (r el ) is its spectral energy distribution ass umed to be 
a Planck function with temperature T gr dBarvainis 1987). The ratio 
between the total luminosity of the hot dust (L™) and that of one 

grain (Lr) gives the number of dust grains ,/Vhd = -pr- was 
obtained by integrating the flux under the Planck function fitted to 
the nuclear spectrum and Lf was obtained by the integration of 
equation[3]for a temperature of T„ = 829 K assuming a = 0.05 //m, 
d ir = 1.4 x 10 24 and y = 1.6 dBarvainis|[l987l ; iKishimoto et al.l 
2007). Finally, the mass of the emitting hot dust is obtained as: 

M w ~ y« 3 AWgr- (4) 

Assuming a graph ite density of p„ r = 2.26gctrT 3 
dGranato & Danes J 19941) . we obtain M H n = (1.4± 0.5) x 10~ 2 M . 
The derived mass is comparable to values obtained for o ther Seyfert 
galaxie s, which are in the range 10~ 6 - 10~ 2 M Q (e.g. Riff el et al.l 
l2009bl l3). 

4.2 Flux distributions 

The flux distributions shown in Fig.[4] are more more extended 
along the PA»135/315° and resemble those in optical lines 
(Bow er et all 1 19951; IStoklasova et alll2009l) . Using HST imaging, 
iBower et all d 19951) observed a narrow "jetlike" feature in the 
[O m]+HyS emission, extending to 174 north-west from the nucleus 
along PA=315°. They concluded that this image is dominated by 
[Om] emission and thus that the high excitation gas is strongly 
concentrated in the jet. To the south-east the [Om] emission is 
much fainter than to the north-w est, similarly to wha t we observe 
in the [Fen] and [Pn] emission. IBower et all d 19951) present also 
a Ho-+[Nn] narrow-band image, which shows the same "jetlike" 
structure, but with a flux distribution broader then in [Om], with 
emission to both sides of the nucleus, similar to our Bry and Pa/? 
flux distributions. The 3.6 cm radio map presents a bipolar jet at 
approximately the same orientation of the emission-line flux distri- 
butions (see Fig.[4]l. By comparing the radio and optical line emis- 
sion, |Brweretal| dl995l) suggest that the low- and high-excitation 
gas present distinct spatial distributions, with the former preferen- 
tially located in the galaxy plane and the latter in the jet region. 
They explain the lower intensity of [O in] to south-west of the nu- 
cleus as being due to obscuration of the jet by the galactic plane. 

A detailed analysis of the line flux distributions reveal some 
differences between lines from distinct ionization levels. In or- 
der to do this comparison, we extracted one-dimensional cuts 
along PA=135° for a pseudo-slit with 0725 width. The resulting 
one-dimensional flux distributions for the [Fen], Pa/?, Bry and 
H 2 /12.12yum emission lines are shown in Figure|7] As observed in 
this figure (see also the two-dimensional maps from Fig.[4]l, the 
[Fen] flux distribution peaks at «0785 north-west of the nucleus, 
approximately coincident with a hot spot observed in the radio 
map, and is fainter to the south-east, in good agreement with the 
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Figure 6. Top: Nuclear and extra-nuclear spectra of Mrkl066. The nuclear spectrum was extracted within a circular aperture with radius C'25, while the 
extra-nuclear spectrum was integrated within an aperture with inner radius ff'25 and outer radius of ff/35. Both spectra were corrected by reddening. The 
extra-nuclear spectrum was multiplied by 1.45, after normalization to the same aperture of the nuclear one, in order to eliminate the stellar absorption in 
the subtracted spectrum. Bottom: Difference between the nuclear and extra-nuclear spectra, together with a fit (solid line) considering the contribution of a 
power-law (dotted line) plus a blackbody function (dashed line). 



[O m] image. The H recombination lines and H2 line emission peak 
at the nucleus. The H lines have a secondary peak at (X'5 south- 
east of the nucleus, which is not observed in the H 2 and [Fe 11] flux 
distributions. The H2 emission presents a secon dary peak to north- 
west, between the nucleus and the radio hot spot. lKnop et all ("2001) 
obtained J and K long-slit spectroscopy of Mrk 1066 with a seeing 
of (y.'7. They present flux measurements for near-IR emission lines 
along PA=45° and PA=135°, which can be directly compared with 
our one-dimensional cuts and flux maps (Figs.|7]and|4j. In general, 
our flux distributions are in good agreement with theirs, although 
much more details are revealed by our measurements. In particu- 
lar, the secondary peaks observed in the H emission lines and the 
H2 nuclear peak are not present in the flux distributions presented 
bv lKnop et aH J200lh . We attribute these differences to our higher 
spatial resolution. 

Our [Fen] and H recombination lines present similar flux 
distributions to t hose of the [Om] and Hq-+[Nii] images of 
iBower et alj Jl995h , respectively. We conclude that the [Fen] emit- 



ting gas extends to high galactic latitudes and is approximately co- 
spatial with the radio jet, while the Bry and Pa/3 emission origi- 
nate mostly from the galactic plane. The H 2 flux is more uniformly 
distributed over the whole IFU field, suggesting also an important 
contribution from gas in the galactic plane, although part of the H2 
emitting gas close to the nucleus may be also outflowing to high 
latitudes, as judged from the highest flux to the north-west than to 
the south-east (Fig. |4j. These conclusions are supported by the ob- 
servation of distinct kinematics for these emission lines, with the 
[Fe 11] originating in a much more disturbed gas than the H and H2 
emitting (Riffel et al., in preparation). This result is also supported 
by near-IR IFU observations of other Seyfert galaxies, which show 
that the H2 and ionized gas have distinct distributions and kine- 
matics, with the former usually restricted t o the galactic plane and 
the lat t er extending to higher latitudes (e.g. lRiffel et al.l2006l . l2008l 
l2009d : lstOTchi-Bergmann et alj|2009l) . 
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Figure 7. One-dimensional cuts from the flux distributions in [Fen], Pa/3, 
Bry and H2/I2. 12/jm emission-lines extracted of a pseudo-slit with (X'25 
width oriented along PA=135°. The fluxes are normalized to the peak value. 



4.3 Diagnostic diagram 

In order to investigate the nature of the circum-nuclear line- 
emitting region we constructed the spatially resolved spectral- 
diagnostic diagram [Feii]/U.25yum/Pa/? vs H2/12.12^m/Bry , 
originally proposed for non-resolved spectra IILarkin et alj 1 19981; 
Rodriguez- Ardila et al. 2004; Rodrfguez-Ardila, Riffel & Pastoriza 



2005a). Typical values for the nucleus of Seyfert galax- 
ies are 0.6<[Fen]/Pa/?<2.0 an d O6<H 2 /Bry<2.0 (e.g. 

iRodriguez-Ardila. Riffel & Pastorizi] l2005al) . The diagram is 
shown in the top panel of Fig.(8] in which black filled circles 
represent typical values for Seyfert galaxies, blue open circles 
for Starbursts and red crosses for LINERs. Most ratios present 
Seyfert-like, although some points are obs erved in the LINER and 
Starburst regions of the diagnostic diagram. Stokl asova et"aH (2009 ) 
analyzed optical IFU data for Mrk 1066 and presented several 
optical diagnostic diagrams. A comparison between the optical and 
the present near-IR diagram shows that they are consistent with 
each other. 

The locations of the distinct line ratios are shown in the bot- 
tom panel of Fig. [8] Seyfert-like line ratios are observed approx- 
imately along the radio jet and ion ization cone (PAx 135/315°, 
iBoweretalJ 19951 ; lNagaretal.il 999 ) suggesting that the line emis- 
sion in these regions is driven by the AGN. The Starburst ratios 
are mainly observed in a region at O'.'S south-east of the nucleus at 
the position of the secondary emission peak of the H recombina- 
tion lines, as described in the Sec. |4.2| Approximately at this same 
position there is also a higher flux in the continuum (see Fig.[2]l. 
Fig.[8]shows also some Starburst ratios at 1" north-west in a region 
approximately coincident with a small enhancement of the H emis- 
sion (see Fig.[7) and with a knot in the optical continuum image 
(see Fig.[2j. We interpret the Starburst ratios as originating in star 
formation regions (SFRs) located at these positions. This conclu- 
sion is supported by the observations of knots in the optical con- 
tinuum image and its absence in the near-IR continuum, as young 
stars contribute mostly to the optical emission, while in the near-IR 
the continuum is dominated by the emission from old stars. At dis- 
tances >Cf!5 perpendicular to the orientation of the ionization cone 
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Figure 8. Top:[Fen]/11.25/jm/Pa/J vj H2/t2.12/jm/Bry line-ratio diagnostic 
diagram. The dashed lines delimit regions with ratios typical of Starbursts 
(blue-open circles), Seyferts (black-filled circles) and LINERs (red crosses). 
Bottom: Spatial position of each point from the diagnostic diagram. 



the H 2 /12.12/im/Bry ratio reaches the highest values, being in the 
region of LINERs in the diagnostic diagram. We interpret this low 
ionization region as being originated by a diffuse radiation field es- 
caping through the walls of the ionization cone, which have enough 
energy to excite the H2. A similar s cenario is observed for the well- 
studie d Seyfert galaxy NGC4151 dKraemer. Schmitt & Crenshaw! 
2008) and attributed to the origi n of the near-IR lines in re gions 
away from its ionization cone Istorchi-Bergmann et alj|2009T) . 



4.4 H 2 excitation 

The molecular hydrogen flux distribution and kinematics in 
the central regions of active galaxies have been the subject 
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Figure 9. Correlation between emission-line fluxes (top left: [Fe n] Al. 2570 jjm, top right: Pa/3, bottom left: [Pn] /11.1886/.im and bottom left: H2 /12.1218/jm) 
and 3.6 cm radio emission. 



of several recent studies. Nevertheless, its excitation mecha- 
nisms is still in debate (e.g. iReunanen. Kotilainen & Prietdl2002l: 
Rodriguez- Ardila et al. 2004; Rodrfguez-Ardila, Riffel & Pastoriza 
2005a| - iDavies et al.ll2005l: iRiffel et alj|2006t |2008j: | Sanchez etal 
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Hicks et al.ll2009l ; IStorchi-Bergmann et al.ll2009h . The H 



2009; 



2 emission 

lines can be excited by two mechanisms: (i) fluorescent excitation 
through absorption of soft-UV photons (912-1108 A) in the 
Lyman and Werner bands - present both in star-form ing regions 
and surrounding AGNs jBlack & van Dishoeckl Il987l) and (ii) 
collisional excitation (usually referred to as thermal processes) 
due to heating of the gas by shoc ks, due to interaction of a ra dio 
jet with the interstellar medium 



Hollenbach & McKee 1989) or 



Malonev. Hollenbach & Tielens 



by X -rays from the central AGN 



We can use the H2/i2.1218/Bry emission-line ratio to in- 
vestigate the origin of the H2 emission. Typical values for Star- 
burst galaxies, where the main heating agent is the UV radia- 



tion, are in the range H 2 /Bry<0.6 teodrfguez-Ardila et al.l l 2004 
iRodrfguez-A rdila. Rif fel & Pastoriza l2005ah . while for Sevferts 
this ratio is larger because of additional thermal excitation by 
shocks or X-rays from the AGN. As observed in Figs. [5] and [8] 
Mrkl066 presents two regions with H 2 /12.1218/Bry < 0.6. The 
first located at ~ 1" north-west from the nucleus and the second 
at «075 south-east. As discussed in Sec . 14. 31 we interpret this re- 
sult as due to an enhancement of Bry due to ionization from young 
stars in SFRs at these locations. In regions away from these SFRs, 
H 2 /12.1218/Bry is larger than 0.6 reaching ~ 2.6 in regions more 
than 077 from the nucleus in the direction perpendicular to the ra- 
dio jet, suggesting additional H 2 emission due to thermal processes. 

Another line ratio which can be used to study the H2 
excitation is H2/12.2477//12.1218. For fluorescent excitation, 
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1 - 0.2 (e.g. Mouril 1994 
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2009). In the case of active galaxies, the fluorescent process 
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Figure 10. Relation between JV U p p = and £ upp = 7", for the H2 emission 
lines for thermal excitation at the nuclear position and at ff.' 5 north-west 
from the nucleus. Ortho transitions are shown as filled circles and para 
transitions as open circles. From left to right the transitions shown are: 1— 
0Q(1), 1-0S(0), 1-0Q(2), 1-0S(1), 1-0Q(4), 1-0 Q(5), 2-1 S(l) and 1- 
0S(2) 



seems not to be important (e. g. Rodriguez-Ardila et al. 20041 : 
iRodriguez-Ardila. Riffel & Pastoriza* l2005al : iRifFel et alj |2006|) . 
As observed in Tab.Q] H 2 ^2.2477/i2.1218~0.1 in all positions, 
indicating that the H2 emission in most of the circum-nuclear 
region of Mrk 1066 is also due to excitation by thermal processes. 

We can also use the observed fluxes for all H 2 emission lines in 
the K band to calculate the thermal exc itation temperature. Follow- 
ing Iwilman, Edge & Juhnsto nj J2005I). we have investigated the 
expression (see also Storchi-Bergmann et all2009T) : 



log 



■ constant ■ 



(5) 



where F, is the flux of the i' h H 2 line, A; is its wavelength, A, is the 
spontaneous emission coefficient, g, is the statistical weight of the 
upper level of the transition, T, is the energy of the level expressed 
as a temperature and T cxc is the excitation temperature. This rela- 
tion is valid for thermal excitation, under the assumption of an or- 
tho: para abundance ratio of 3:1. In Fig.[K)]we present the observed 
values for _/V U pp = jj- (plus an arbitrary constant) vs E upp = T, 
for the nuclear spectrum and for an extra-nuclear one at 075 north- 
west from the nucleus (position A in Fig.[T]l. The resulting fit of the 
above relation is shown in Fig.[l0]as a continuum line and resulted 
in an excitation temperature of T exc = 2131 ± 40 K for the nucleus 
and r cxc = 2097 ± 40 K for the position A. The fact that the or- 
tho lines (1-0 S(l), 2-1 S(l), 1-0 Q(l), etc) give the same result as 
the para lines (1-0 S(0), Q(2), Q(4), etc) confirm that the H 2 is in 
thermal equilibrium at T cxc and that the ortho to para ratio is 3 as 
assumed. This is also what is expected for thermal equilibrium. 

The heating of the H 2 can be provided by nuclear X-rays 
and/or by shocks due to the interaction of the radio jet with the in- 
terstellar medium. Nevertheless distinguishing between both mech- 
anisms is usually a hard task. Correlations between the radio or X- 
ray emission with the H 2 emission c ould be useful to inv estigate 
the importance of these mechanisms. lOuillen et alj d 1 999t> . using 
HST, imaged a sample of 10 Seyfert galaxies in H 2 and looked for 



correlations with radio 6 cm and hard X-ray flux. They found no 
correlation with X-rays and a weak correlation with radio 6-cm, 
suggesting that no single mechanism is likely to be responsible for 
the H 2 excitation in Seyfert galaxies. 

In Fig.[9] we present plots for the [Fen] /ll.2570/jm (top 
left panel), [Pn] /il.l886yum (bottom left), Pa/3 (top right) and 
H 2 /12.1218/mi (bottom right) emission-line fluxes versus the radio 
flux at 3.6cm for each spatial pixel. As observed in this figure the 
flux distributions correlate very well with the radio emission. We 
used the routine r.correlate in idi0 programing language to obtain 
the Sepearman correlation coefficient (R). The best correlation is 
found for [Fen] with R = 0.77, followed by Pa/3 with R = 0.70, 
[P 11] with R = 0.64 and H 2 with R = 0.62. These correlation can be 
understood in two ways: (i) the interaction between the radio jet and 
the ISM produces the heating necessary to collisionally excitation 
the emitting gas and/or (ii) the radio jet only compresses the gas 
enhancing its emission, excited by other mechanism. The fact that 
Pa/8 emission shows a better correlation with the radio than the H 2 
emission suggests that the radio jet only compresses the gas, since 
the H recombination lines are not sensible to collisional excitation. 
Nevertheless, shocks due to the radio jet cannot be completely dis- 
carded since the ionized and molecular hydrogen emission could 
arise from distin ct locations along the lin e -of-sig ht. 

Following IRodriguez-Ardila et al.l d2004f> and iRiffel et al.1 
d2008h we evaluate the emergent flux for the H 2 /12.1218/um of a 
gas cloud illuminated by a source of har d X-rays using the model of 
iMalonev. Hollenbach & Tielensl ( fl996T) . IShu et alj d2007l) obtained 
a 2-10 keV X-ray flux of F x ~ 2.3 x 10~ 13 ergs -1 cm -2 and a attenu- 
ating column density of N H > 10 24 cnT 2 for Mrk 1066 using Chan- 
dra observations. The N H value is measured along the line-of-sight 
and according to the unified model N H may be smaller in other di- 
rections, since Mrk 1066 has a Seyfert 2 nucleus and thus the dusty 
torus is seen ed ge-on, enhancing the column density only along 
the line-of-sight (Antonucci 1993). The column density in the H 2 
emittin g region can be estimated usin g Nh ~ 5.2 x 10 21 E(B — V) 
[cm- 2 ] Ishull & van Steenberdl 1985b . Using the E(B - V) values 
from Fig. [5] we obtain N H =3 1 21 cm -2 , which is in good agre e- 
ment with the value obtained bv lGuainazzi. Matt & Perolal d2005l) - 
N H = 1.21 x 10 2 ' cirr 2 - using XMM-Newton X-ray observations. 
The lower value of N H obtained with the XMM-Newton relative to 
the one with Chandra data may be due to the fact that former has a 
larger aperture than latter, diluting the effect of the nuclear dust. Us- 
ing N H x 10 21 cnr 2 as the appropriate value for the emitting gas, 
the emergent fluxes calculated from the X-ray excitation model can 
account for the H 2 emission in most locations of the observed field 
of Mrk 1066. 

The line H 2 2-1 S(3) (/12.0649yum) can also be used to investi- 
gate the H 2 excitation. In the case of X-ray i rradiated gas this line 
is expected to be absent (Daviesetal. 2005 ). Our spectral range 
does n ot include this line, but Ri ffel. Rodriguez-Ardila & Pastoriza 
(2006) obtained a near-IR nuclear spectrum for an aperture of 
l'.'6x076 in which the H 2 2-1 S(3) emission line is not detected, 
supporting the above conclusion that X-ray heating contribute to 
the observed H 2 emission. 

From the discussion above, we conclude that heating by X- 
rays from the central AGN may be the dominant excitation mecha- 
nism of the H 2 . This conclusion is in good agreement with those ob- 
tained for other active galaxies (e.g.lStorchi-Bergmann et alj 19991 : 
IRiffel et al.ll2008l : IStorchi-Bergmann et al.ll2009l) Nevertheless, as 
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discussed above, some contribution of the radio jet to the H 2 excita- 
tion cannot be discarded. In fact, the interaction of the radio jet with 
the circum-nuclear gas seems to be important in the excitation o f 
the H 2 in some galaxies, such as ES0 428-G14 (Riffe l et al]2006l) . 



4.5 [Fen] excitation 

The [Fen] /il.2570yum/Pa/3 line ratio is controlled by the ratio be- 
tween the volumes of partially to fully ionized regions, as the 
[Fe n] emission is excited in partially ionized gas regions. In AGNs , 
such regions can be created by X-ray (e.g. ISimpson et all 1 1996T) 
and/or shock (e.g. lForbes & Ward 1993) heating of the gas. 

The [Fe n] /ll.2570/im/Pa/? line ratio can be used to in- 
vestigate the excitation mechanism of the [Fen] emission. 
For Starburst galaxies, [Fen]/Pa/?< 0.6 and for supernovae 
for which shocks are the main excitation mechanism, this 
ratio has values larger than 2 I Rodrfguez-Ardila et alj 1 20041 ; 
Rodri guez- Ardila. Riffel & Pastorizall2005af) . Thus this ratio can 
be used to measure the relative contribution of photoionization 
and shocks, with values close to 0.6 meaning that photoionization 
dominates, while values close to 2 that shock excitation is dom- 
inant. As observed in Figs.[5] and [8] for Mrk 1066 typical values 
for this ratio are 0.6-1.5, although in the regions co-spatial with 
the SFRs [Fe ii]/Pa/3< 0.6 due to an enhancement of the Pa/3 emis- 
sion. The highest values are observed in locations close to the edge 
of the radio structure suggesting that the excitation due to shocks 
by the radio jet is important in these regions. The stronger cor- 
relation between the [Fe n] and the radio emission than those ob- 
served for the other emission lines (see Fig.[9} support a higher 
contribution of the radio jet for the excitation of [Fen] than for 
H2. Radio shocks have been appointed as the main [Fen] exci- 
tation mechanism for some Seyfert galaxies such as ES0428- 
G14 teiffel et alj 2006) and being not ne gligible in others, as in 
NGC4151 Jstorchi-Bergmann et al]2009l) . 

Another tracer of the origin of the [Fen] emission is the 
line ratio [Fen]/1 1.2570 pm/[Pn]/l 1.8861 pm. These two lines have 
similar excitation temperatures, and their parent ions have similar 
ionization po tentials and radiat ive recom bination coefficients. As 
pointe d out in lOliva et alj J2001I) - see also lStorchi-Bergmann et al.1 
(l2009h - values larger than 2 indicate that shocks have passed 
through the gas destroying the dust grains, releasing the Fe and 
thus enhancing its observed abundance. For supernova remnants 
where shocks are the do minant excitation m echanism [Fe n]/[P 11] is 
typically higher than 20 JOliva et alj|200lh . As observed in Fig.[5] 
Mrk 1066 presents [Fen]/[Pn]~3 at most locations, but in some re- 
gions close to the borders of the radio structure [Fe n]/[P 11] reaches 
values up to 9.5, indicating that shocks due to the radio jet are im- 
portant in these locations in agreement with the highest values ob- 
tained for the [Feii]/Pa/3 at the same locations. 

From the discussion above, we conclude that X-ray heat- 
ing may be the dominant excitation mechanism of the [Fen] in 
Mrk 1066, exceptions are the regions tt'85 north-west of the nu- 
cleus where there is a radio hot spot and in the borders of the radio 
structure, where shocks seem to play a more important role. In any 
case, the contribution of shocks due to the radio jet seems to be 
more important for the [Fe 11] emission than for the H 2 emission, as 
indicated by the better correlation observed for with the radio emis- 
sion and the [Fen] than the one observed for the H2. This result is 
in good agreement with the ones obtained for other galaxies (e.g. 
Rodriguez- Ardila et al. 2004; Rodriguez- Ardila, Riffel & Pastoriza 
2005a; Rif fel et alj2006l ; IStorchi-Bergmann et aljl999ll2009h . 



4.6 Mass of ionized and molecular gas 

Following iRiffeletalJ J2008h and IStorchi-Bergmann etal] J2009h 
we estimate the mass of ionized gas in the inner 700 x 700 pc 2 by 

M HU « 3 x 10" ( F f 2 | f-^f [M ], (6) 
\ erg s 1 cm 1 J \ Mpc / 

where F Bly is the integrated flux for the Bry emission line and d is 
the distance to Mrk 1066. We have assumed an electron tempera- 
ture T = 10 4 K and electron density N e = 100 cnT 3 . 
The mass of molecular gas can be obtained as 



M„. ~ 5.0776 x 10 1: 



H 2 .12.1218 



erg s -1 cm 2 / \ Mpc 



(- 



[Mo], 



(7) 



where F Hl n.ui» is the integrated flux for the H2 /I2.1218yum emis- 
sion line and we have used the vibrational temperature T=2000 K 
we have obtained in Sec. 14.41 

Integrating over the whole IFU field we obtain F^ Ty w 2.5 x 
10~' 4 erg s -1 cm 2 and -F tf2 ^ 2 .i2i8 ~ 2.8 x 10~ 14 erg s -1 cnT 2 , result- 
ing in Mmi ~ 1.7 x 10 7 M o and M Hj a 3.3 x 10 3 M o . The mass 
of molecular gas is 10 4 times smalle r than of ionized gas but, as 
discussed in I S torchi-B ergmann et al. ( 2009), this H 2 mass repre- 
sents only that of hot gas emitting in the near-IR, heated by X-rays 
and maybe also by shocks from the central AGN. The total mass 
in molecular gas (including the cold gas which do es not emit in the 
near-IR) is usually 10 5 -10 7 times that in hot H 2 jDale et alj|2005h 
suggesting that the total molecular gas mass is at least 10 8 M . 



5 CONCLUSIONS 

We used integral field J- and KpBand spectroscopy of the inner 
700 x 700 pc 2 of the Seyfert galaxy Mrk 1066, obtained with the 
Gemini NIFS at spatial resolution ~35pc and spectral resolution 
~3.3 A, to map the near-IR continuum and emission-line flux dis- 
tributions, as well as the properties of the nuclear spectrum. Our 
main conclusions are: 

• The nucleus contains an unresolved infrared source whose 
continuum is well reproduced by the emission of a dust structure 
with temperature ~830 K and mass ~ 1.4x 10~ 2 M Q which we iden- 
tify as the dusty torus of the Unified Model; 

• Emission-line fluxes (except for the coronal lines) are elon- 
gated in PA= 135/315° in agreement with previous optical [Om] 
imaging and also following the radio jet extending beyond the IFU 
field. Except for the H lines, the emission is stronger to the north- 
west than to the south-EAST in association with the radio hot spot 
implying that at least part of the emitting gas is co-spatial with the 
radio outflow. The H emission is stronger to the south-east, where 
we find a large region of star-formation. 

• There is a strong correlation between the radio emission and 
the highest emission-line fluxes suggesting that the radio jet plays 
an important role for these intensity levels. At lower line flux values 
this correlation disappears indicating a contribution from gas in the 
galactic plane. 

• The coronal line emission in [Ca vm] /12.3220yum and 
[Six] A 1.2524 yum are unresolved by our observations indicating 
that they originate within a radius of 18 pc from the nucleus. 

• The [Fenj/Pa/? vs H 2 /Bry line-ratio diagnostic diagram 
presents typical values for Seyfert nuclei along the ionization cone 
indicating that the line emission is powered by the AGN in these 
regions. Some radiation escapes through the walls of the ionization 
cone and powers a low ionization emission-line region in locations 
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perpendicular to the cone and radio jet. At the positions Of.' 5 south- 
east and at 1" north-west of the nucleus Starburst-like values are 
observed due to additional emission of the H lines from star form- 
ing regions. 

• The reddening map obtained via the Pa/3/Bry line ratio 
present a S-shaped structure with high values oriented along the 
PA«135/315°, reaching up to E(B - V) = 1.7. Lower values are 
observed in regions away from this structure with typical values of 
E(B - V) = 0.5. The reddening in the continuum is smaller with 
typical values of E(B — V) = 0.9 along the S-shaped structure. 

• The H 2 gas has an excitation temperature r exc « 2100K and 
its emission is mainly due to X-ray excitation from the central 
AGN, which is also the main excitation mechanism of the [Fen] 
emission. Emission due to shocks produced by the radio jet may 
contribute a small fraction to the emission of both lines, with a 
larger contribution for the [Fe n] line. 

• The mass of ionized gas in the inner 700x700 pc 2 of Mrk 1066 
is M HU a 1.7 x 10 7 Mq, while that of hot molecular gas is M Hl x 
3.3 x 10 3 M . 
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